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ABSTRACT 

In this letter we introduce a new method to determine the radiative neutrino 
decay rate in the presence of a medium. Our approach is based on the gene¬ 
ralisation of the optical theorem at finite temperature and density. Differently 
from previous works on this subject, our method allows to account for dispersive 
and dissipative electromagnetic properties of the medium. Some inconsistencies 
that are present in the literature are pointed-out and corrected here. We shortly 
discuss the relevance of our results for neutrino evolution in the early universe. 

1 Introduction 

Most of neutrinos in the universe are expected to be produced in the presence of 
environments characterised by very high temperatures and densities. In such media 
neutrinos behave quite differently from the way they do in vacuum jl]. The coherent 
interaction of neutrinos with the particles belonging to the medium can amplify or 
induce effects that are otherwise very feeble or absent in the vacuum. This can be 
crucial in order to disentangle possible anomalous properties of the neutrinos that 
may testify for new physics present beyond the standard model. MSW |2j neutrino 
resonant oscillations give a well know example of such a kind of effects. 

Neutrino electromagnetic properties are also expected to be significantly affected 
by the presence of media. In fact, due to their coherent interaction with charged 
leptons and nucleons in the medium, neutrinos acquire an effective coupling to the 
electromagnetic field. This effect was investigated by several authors who showed 
how even massless neutrinos passing through the matter acquire an effective charge 
PI and, in a charge asymmetric medium, also an effective magnetic dipole moment 

1 New address since 1st October 1998: Dipartimento di Fisica, Universita di Padova, Via F. 
Marzolo 8, 1-35131 Padova, Italy 
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[Q. Several intriguing consequences of medium induced neutrino electromagnetic 
couplings have been studied in the literature, e.g. plasmon decay in vv pairs || and 
Cherenkov emission of neutrinos |J. 

In this letter we focus on the radiative decay of a heavy neutrino into a lighter one. 
Our underling particle physics model consists of a minimal extension of the standard 
model which allows for non-vanishing neutrino masses and mixing. In the ambit of 
this model, neutrino radiative decay is extremely suppressed in vacuum because of the 
Glashow-Iliopoulos-Maiani (GIM) cancellation fTj. However, if T < m T , this is not 
the case in a heat-bath due to the different thermal populations of the three charged 
leptons families. As a consequence, the radiative neutrino decay rate in the presence 
of hot media exceeds the decay rate in vacuum by many order of magnitudes. This 
was first showed by D’Olivo, Nieves and Pal (DNP) || (see also @), who computed 
the decay rate both in the case of a non-relativistic plasma (NR) (T -C m e ) and of 
a ultra-relativist plasma (UR). Recently, DNP’s work has been extended by Nieves 
and Pal |I0| (NP) who showed that the decay rate may be further enhanced due to 
the Bose-Einstein stimulation in the production of low-energy photons (a; « T). 

Although we agree with DNP and NP concerning the huge neutrino decay rate 
amplification occuring in a medium we do not , however, about their method and some 
of the results obtained in their works. Our main disagreement concerns the properties 
of the photon produced by the neutrino decay. This photon was assumed to be on the 
light-cone (A' 2 = 0) by DNP and NP. Such an assumption is equivalent to disregard 
any dispersive and dissipative property of photons in a thermal bath. It is well known, 
however, that in plasma photons behaves more like collective excitations (usually 
dubbed plasmons ) than like free particles. While neglecting collective properties of 
photons may be a reasonable approximation in the high frequency regime this is 
however inconsistent when they are low energy or, so-called, soft. We will show that 
to account for the correct plasmon dispersion relation is indeed mandatory in order to 
avoid an unphysical behaviour of the decay rate as a function of the neutrino velocity. 

Photon dissipative properties need also to be taken into account. In fact, in a heat 
bath low-energy photons are continuously absorbed and re-emitted mainly by Landau¬ 
damping on the heat-bath free-charges. For the first time we will show that Landau¬ 
damping gives rise to an additional, and in many cases dominant, contribution to 
the decay rate. Since the photon asymptotic state is a not a well defined quantity 
in a medium, we recourse here to a treatment that is based on finite temperature 
field theory and on the generalisation of the optical theorem to such a framework 
[ |TT] . Within this framework the basic quantities are the Green functions which are 
well defined quantities even at finite temperature and density. As a bonus, one get 
that using such approach the role of Landau-damping comes-out naturally from the 
thermally-corrected photon progator |12|. A similar method has been previously 


applied by other authors to determine the chirality-flip rate of Dirac neutrinos in a 
degenerate plasma ||T3|, |T4| and in the early universe (15| . 
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Figure 1: Feynman diagram contributing to the heavy neutrino self-energy E. 


In this letter we only consider the case of an UR non-degenerate electron-positron 
plasma. As we will shortly discuss in our last section, this case is relevant for the 
study of neutrino evolution in the early universe. The cases in which the medium is 
degenerate and/or NR will be discussed elsewhere |L6|. Our general method will be 


presented in section 2. Section 3 contains the computation of the decay rate in the 
case in which the plasmon is on-shell. The generalisation of this result to the case 
with off-shell photons is discussed in Sec. 4. Sec. 5 contains a discussion about the 
applications of our results to the early universe physics. 

2 The general method 


It was first showed by Weldon |1T[ that in a heat-bath the decay rate of a fermionic 
particle species (having the energy E ), and the rate of the inverse process Fj, are 
related to the imaginary part of the fermion self-energy E by the following relation 


Im {w ( p)E(p)w(p)} = -E(Ti + r d ) 


( 1 ) 


This result is the generalisation of the optical theorem for finite temperature and 
density. It is understood that E includes corrections due to the effects of the medium 
on the held propagators. 

In our case, the interesting contribution to the heavy neutrino self-energy comes 
from the Feynman diagram represented in Fig. 1. In that figure zy and Uj are, 
respectively, the heavy and the light neutrino mass eigenstates participating to the 
decay and inverse-decay processes. Working in the real-time-formalism (RTF) we 


determine the imaginary part of E by properly cutting [17j the diagram represented 
in Fig. 1. We get 

e(p 0 ) f d 4 K 


Im E (p) = 


2 sin 2 (j) p J ( 27 t ) 4 


r“*(A') S-(p-K)D*pK) + S*(p-K)D-JK) T„(K) 


/IV v 


( 2 ) 
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where 


1 e ^bol/2 

2 rp e P\Po\ +1 


( 3 ) 


and the expression for the effective vertex Y^(K) will be given below. Here K = 
(fc 0 , k) is the photon 4-momentum in the medium rest-frame and e(po) is the sign 
function. We observe that a positive k 0 accounts for the process with emission of a 
photon (energy given to the medium) whereas a negative k o takes care of the inverse 
process (energy lost by the medium). We are then able to disentangle the different 
contributions to Y^ and Ik from (jl) by properly choosing the k 0 integration interval. 
For the sake of simplicity we assume the light neutrino to be massless. In the RTF 
the off-diagonals component of its propagator are 


S ± (p - K) = -27 ri (p - K) (0(±(p o - k 0 ) - n F (\p 0 - fc 0 |)) S ((p - K ) 2 ) (4) 


where n F {B) is the Fermi-Dirac (Bose-Einstein) distribution function. In the rest 
frame of the medium (and in the Landau gauge) the photon propagator can be de¬ 
composed into 

D%(K) = D$P,„(K) + DtQ„„(K) . (5) 

The transverse and longitudinal projector operators are defined by 


PAK) 



kikj 

~k 2 " 


3kj 


Q^{K) = e^(K)ekK) 


( 6 ) 


where e (U (iF) = (k, k 0 k / k) / \JK 2 is the polarization versor. The coefficients L are 
given by p 


D fl = 2 (6(±k 0 - n B (|fc 0 |)) 


nS (. K 2 - Re Yi TL {K)) 5 Kr (lmU T . L (K)) 


_M rjjjq _ 

(K 2 - ReU TtL (K)f + (Imn t, L ( K )) 2] 


(7) 


where 5 F r is a Kroneker function. These expressions include the one-loop thermal cor¬ 
rections to the photon propagator. It is known JTD|] that for soft photons (ko, k < eT) 
these corrections cannot be neglected and that a meaningful perturbative expansion 
can be obtained performing a proper resummation of the hard-thermal-loops con¬ 
tributing to the photon polarization tensor. 

Equation ([?]) has been written in such a way in order to distinguish an on-shell 
part, corresponding to the propagation of the transversal and longitudinal modes of 
a plasmon, and an off-shell part associated to Landau-damping. The fluctuation- 
dissipation theorem shed light on the physical nature of off-shell plasmons: they are 
thermal fluctuations of the electromagnetic field [IT3]. As suggested by the Breit- 


Wigner form of the off-shell part of Df L , we can think of electromagnetic thermal 
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fluctuations as resonances which are continuously emitted and re-absorbed from the 
free-charges in the plasma by Cherenkov and inverse-Cherenkov emission. 

The general expression of the polarization tensor in a isotropic plasma is 


IWAT) = n T (K)P^(K) + n l (k)q^(k) 

where, for an UR plasma, we have [ 20[ 

k$ -\- k 


( 8 ) 


Re II T (iC) = 


3up 


M + (1 _M)h hl 

k 2 [ k 2> 2k 


k n — k 


lie 11,.(A') = 3u>p(l — f|) 


1 - h in 

2k 


k 0 + k 


kn — k 


(9) 


and 


lmU T (K) = 


3 u>p . fcfl, k 0 


2 

im n l (k) = 3wp7r(i - h) h _ fcg) . 


( 10 ) 


In the UR limit, the plasma frequency is given by Up = eT/ 3. 

Plasmon dispersion relations are determined by solving K 2 = Re I\r.riK). Ana¬ 
lytic forms of the solutions can be obtained in some suitable limits [12[ . In the k —> 0 
limit transversal and longitudinal dispersion relations are respectively 

6 
-I 

5 
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k 2 = u% + -k 2 transversal , 


( 11 ) 


k 2 = Up + - k 2 longitudinal , 


whereas in the hard limit (fc 0 , k^> up) one gets 
kl = \oj 2 p + k 2 


transversal , 


k 0 = k 2 + 4 k 2 e 2k2 / 3uJ p longitudinal . 


( 12 ) 


3 The on-shell plasmon case 

We now come into some more details of our computation. As we discussed in the 
introduction, in a medium the neutrino-photon coupling is mediated by the weak 
interaction of the neutrino to the free charges in the plasma. We the sake of simplicity 
we only consider here a non-degenerate electron-positron plasma at temperatures 
below the muon mass. From [[3]. |2T| we know that the effective medium induced 
neutrino electromagnetic vertex is given by 

r p(K) = U ie —Hp„(K)YL (13) 

e 
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where U is the lepton mixing matrix and L = \{1— 75 ) 


By substituting (p~3|) , (|4j) and the on-shell part of (0) in (^), we find that the decay 
rate of the heavy neutrino into the lighter one and a transversal plasmon is 


-pon _ 

1 T — 


1 


2 sin 2(p p 


Gl ur*rr ,2 f^K 

4 Ee 2 ' ie jel J ( 2 vr ) 4 
x [e(p 0 - k 0 )(6(-k 0 ) + n F (\p 0 - fc 0 |)) + e(p 0 - k 0 )(6(k 0 ) + n B (\k 0 \))} 


(ReU T (Ii )) 2 S (K 2 - Ren T (A')) 5 ((p - IiJ 


xP Q/3 Tr {(p — m) 7 ^(p — Ji) 7 q l| 


(14) 

where m is the heavy neutrino mass. In the case the plasmon is longitudinal the rate 
can be obtained by replacing in (|T3|) the subscript T with L and the projector P a g 
with Q a p. In the square brackets in the right side of (|T|) we can distinguish a term 
accounting for the Bose-Einstein stimulated photon emission and a Pauli-bloeking 
term for the production of the light neutrino . Since generally the latter is much 
smaller than the former we disregard the Pauli-blocking term in the following. We 

also approximate the factor-— with the unity. 

sm 2 <p p 

Using 5 ((p — K ) 2 ) and 5 (. K 2 — Relh^/t)) to suppress respectively the angular 

, , , f , d 4 K 

and the temporal parts of the integration we get 


(2tt) 


r° T n = 

32ne 2 


9G2p I u; e a 


je\ 


2 UJpf T {v) 


(15) 


Here v is the velocity of the heavy neutrino with respect to the medium rest frame. 
Since it is in general impossible to find an analytical form of the plasmon dispersion 
relation, it is convenient to distinguish between a soft and a hard contribution to 
In the case the plasmon is transversal we have 


/r f » ^ ^ 

and 


Ml 0 = fr» + 

f kmaX dk{l+n B (l)) \e 2 - -(1 + m 2 ) - i 
9 m 2 v Jk%* 1 1 l 2 l ’ k 2 L 


E 


2(1 +^ 2 ) 


fp d (v) ~ 


iuhard 
I ™max 


—, id+„»(*))Ri-I « ( i-D 

m‘v Je~:' k 2 44 


(16) 

61(1—fc) 

(17) 

(18) 


where E is the energy of the heavy neutrino and tilted quantities have been normalised 
to (Up so that fx(v) is a dimensionless function. The integration limits are 



^hard 

W 2 V 1 - V 2 



(19) 
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Figure 2: In this figure the dashed, dotted and continuous lines represents respectively 
the functions /£ ard (u), /|? ft (r>) and the sum of the two. Here we fixed fh = 5. 


in the hard limit, and 

k s ° f * = rh\ 


1 1 + v 2 

1 — V 2 


- 1 


/.soft _ 
™min 


m\ 


'l-v 2 

1 + V 2 


-1 


( 20 ) 


in the soft one. For reason of available space, only a zero order expansion in k of 
Ren^/i) has been reported in the expressions (|1^) and ([T7|). However, a more 
complete second order expansion has been used to get Fig. 2. Also, in our numerical 
integrations some more suitable exponential cut-off functions have been used instead 
of the step functions appearing in the integrals in (|1^) and (|T7|). From Fig. 2 we 
see that in spite of the large enhancement of the Bose-Einstein factor ns(|A)o|) = 
( e dh'o| _ i)- 1 ~ T/ko 1 that one gains in the soft limit, the soft contribution to 
the total decay rate is subdominant with respect to the hard one. This is due to 
the smaller available phase space when the plasmon is soft. Not surprisily, such a 
situation becomes even more pronounced for higher values of the decaying neutrino 

mass. The maximum value of /| oft 

7 soft _ n 
™min u ‘ 

It is worthwhile to observe that, a part for the presence of the infra-red cut-off 


,i i i • m 2 ~ u p 

i T is reached at the velocity u* = —-- T 

m 2 + ujp 


where 


function (and of an extra factor of 2), our expression (|18|) coincides with the result 
found by NP. In our opinion, the introduction of the cut-off is essential in order to 
avoid the incorrect use of the hard plasmon dispersion relation ( K 2 ~ 0) in the soft 
limit. In fact, in the limit in which the decaying neutrino is ultra-relativistic (v —> 1) 
it is evident from (IHf ) that —> 0. As we have seen above, our approach solves 

such a problem by properly separating the hard and the soft contributions to 
To confirm us of the validity of our method, we get a physically correct behaviour of 
the decay rate that drops to zero when v —> 1 (see Fig. 2). Such a result is actually 
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convincing since no rest frame is available for the neutrino to decay in when it travels 
at the speed of light. A non-vanishing rate was instead found by NP in the limit 


v —> 1. 

In the case the plasmon is longitudinal only the soft term contributes to the total 
decay rate. In fact, as we can see from the second of the equations (0 A' 2 , hence 
also Rell^(A'), vanishes in the hard limit. Due to angular momentum conservation, 
Tl(K) is suppressed by a factor (c op/m) 2 with respect to the soft contribution to 
r-r(A'). For this reason its contribution to the total decay rate is subdominant. For 


reason of available space the expression of Tl(K) will be presented elsewhere [|16 


4 The off-shell plasmon case 

In this section we discuss the contribution of thermal fluctuations of the electromag¬ 
netic field to the neutrino decay. As we discussed in the Sec. 2 the effect of these 
fluctuations is accounted by the second term on the right side of (^). Differently 
from the case considered in the previous section, plasmons associated to thermal 
fluctuations are not on-shell hence they do not obey any dispersion relation. As a 
consequence, in these case we have to deal with a double integration both over k and 
k$. Indeed, the expressions for the decay rate look in this case 


r of! _ 9Gf I 
T - L 32ne 2 


KU je \W P g T , L (v) 


( 21 ) 


The expression of the kinematical functions gr,L(v) are 


9t{v) ~ 2 T- 


m 2 v 


dL 


nsoft 
' ""max 


o L 


soft 


dk 

F 


E 


(A' 2 + m 2 ) (Ek 0 - \{K 2 + m 2 )) 


k 2 


x K 2 A T (k 0 ,~k)9(k 2 -k 2 ) 


( 22 ) 


and 


9l(v) ^ 


where 


T 1 


2 m 2 v Jo 


dk 


dko / — 

ifc soft A; 4 

m i n. ' v 


2E 



x K 4 A L (k 0 , k)6(kr 


k 2 0 ) 


~ ~ (ReU T , L (K)) 2 + (lmU T}L (K)) 2 

A-T,L\K 0 ; = 


(k 2 


m 


(K 2 - Ren T , t (A')) 2 + !lmn r ,,!/rif 


(23) 


(24) 


Since Landau-damping is effective only for soft photons we can neglect here any 
contribution coming from the hard part of the photon spectrum. This allowed us to 
perform the substitution (1 + ns(fco)) —»> T/ko in (|22| , |23|) . The functions gr,L(v ) have 

















Figure 3: In this figure the dashed and the continuous lines represents respectively 
the functions gr(v), 9 l(v). Again, we fixed here m = 5. 


been computed numerically and their plots are reported in Fig. 3. From this figure 
we see that like for the soft contribution to reach their common maximum 

value when v = v*. This is not unexpected since the functions At,l(K) take their 
maximum values when K 2 approaches Re YItj (K). 

We have, however, a crucial difference between the off-shell and the on-shell con¬ 
tributions to the total decay rate. Whereas radiative neutrino decay is kinematically 

forbidden if K 2 > m 2 this is clearly not the case whenever the photon is off-shell. 

3 

Note that even in the hard limit, where K 2 ~ —Up (see (|12])), the on-shell decay 

3 

cannot take place unless m 2 > -Wp. This was not noted in || ||, |I0| . We have to 

mention, however, that although T°^ remains different from zero it is typically very 
small when m -C oop. 


5 Discussion 


The results of our previous sections can be directly applied to study neutrino evolution 
in the early universe where, for T 1 MeV, the electron-positron plasma is UR and 
non-degenerate. For definitess we consider the radiative decay of a tau-neutrino with 
a mass in the MeV range into a massless v e . In order to determine if radiative decay 
can give rise to any relevant depletion of tau-neutrinos we have to compare the decay 
rate with the expansion rate of the universe. Using our previous results we find that 
the decay rate is given by 


T ~ 2 x 10 


-5 


T 


1 MeV 


sin 2 2d | 


f(v)+g(v) 

10 


-i 


(25) 
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where 9 is the vacuum mixing angle between the v e and the v T and f(v) and g(v) are 
the sum of the transversal and longitudinal kinematical functions computed in the 
previous sections. As we showed above f{y ) + g{v) is typically of the order of few 
tens if m ~ T. At the big-bang nucleosynthesis (BBN) time the Hubble rate is given 
by 


H ~ 1 


T \ 2 
1 MeV / 


-i 


(26) 


It is easy to verify that T can exceed H only for temperatures that are above (though 
not too much) the range at which BBN take place (0.1 < Tbbn 10 MeV). Fur¬ 
thermore, even if T > H at such high temperatures, one should not forget that no 
neutrino depletion could take place since, for T > 2m, standard weak processes are 
always able to re-produce tau-neutrinos. 

A different scenario might be realised in the case one considers a larger exten¬ 
sion of the standard model which allows for the presence of a neutrino transitional 
magnetic moment. A magnetic moment induced decay can become dominant with 
respect to that induced by the medium when the plasma in NR and neutrinos are 
out of thermal equilibrium. Although Landau-damping gives an extra contribution 
to the decay also in this case, we do not expect that this effect could give rise to 
a significant neutrino depletion after BBN. We observe, however, that astrophysical 
and cosmological constraints on the neutrino transitional magnetic moments, that 
have been previously derived disregarding dispersive and dissipative properties of the 
medium, should be reconsidered in the light of the results obtained in this work. This, 
as other related issues, will be discussed in a forthcoming paper 
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